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Matter-neutrino resonances (MNR) can occur in environments where the flux of electron an¬ 
tineutrinos is greater than the flux of electron neutrinos. These resonances may result in dramatic 
neutrino flavor transformation. Compact object merger disks are an example of an environment 
where electron antineutrinos outnumber neutrinos. We study MNR resonances in several such disk 
configurations and find two qualitatively different types of matter-neutrino resonances: a standard 
MNR and a symmetric MNR. We examine the transformation that occurs in each type of resonance 
and explore the consequences for nucleosynthesis. 
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I. INTRODUCTION 

Neutrinos shape the physical phenomena surrounding compact object mergers, from the dynamics of the disk or 
hypermassive-neutron star itself [M, to the energetic jets, e.g. @ that mw from them. Neutrinos also play an 
important role in the nucleosynthesis that takes place in and around disks [6l-l9|. For example, the wind outflows 
[id . [Tljl above disks can be home to nucleosynthesis, including perhaps the r-process, depending on neutrino flavor 
composition 0,1 EMI- The neutrino flavor composition above the neutrino trapping surface depends not only 
on thermodynamics in the trapped regions, but also the oscillation of neutrinos as they leave the disk. The high 
neutrino density coupled with high matter density provide an environment where several kinds of oscillation may take 
place. Neutrinos emitted from mergers can undergo the same types of transformations that neutrinos from supernovae 
do [la, as well as oscillations not previously seen elsewhere (except for in collapsars [3), called Matter-Neutrino 
Resonance (MNR) transitions M- The MNR takes place when the matter potential and the neutrino self interaction 
potential are the same size and have opposite signs. 

Much of the previous work on neutrino transformation in high neutrino density environments has considered matter 
and self interaction potentials of the same sign. For example, for the case of core collapse supernovae neutrinos, it has 
been pointed out that several types of transformations can occur in such systems which have a slight electron neutrino 
excess, start at high neutrino density and end at low neutrino density; for recent work see e.g. |20l - l27| |. At very high 
densities of neutrinos, synchronized neutrino oscillations, e.g. [^ . change neutrino flavor on a small scale. At very 
low densities of neutrinos, the neutrino self interaction potential is unimportant and Mikheyev Smirnov Wolfenstein 
(MSW) oscillations can take place when the scale of the matter potential is the same as the vacuum scale. In 
between these extreme regimes of synchronized and MSW oscillations, large scale flavor transformation can take place 
when the neutrino self interaction potential approaches the vacuum oscillation scale and the neutrinos and 

antineutrinos enter the transition region nearly in flavor eigenstates, e.g. [29l - l3^ . During the transition, the neutrinos 
and antineutrinos are said to be “locked” as their survival probabilities mirror each other and the phenomenon is 
referred to as “bipolar” or “nutation” oscillation. 

Unique to settings where the neutrino interaction potential and the matter potential have opposite sign, another 
oscillation phenomenon can be possible [mil- In [I^ it was shown that collapsar-type disks may be home to MNRs 
and that the MNR may result in flavor transformation that alters nucleosynthesis. The understanding of the MNR 
was expanded in M , where the standard MNR was explored in detail. A MNR begins when the scale of the matter 
potential and the neutrino self interaction potential are the same and can cancel. A MNR transition can cause a 
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dramatic change in the flavor of neutrinos. During the transformation, the neutrino self interaction potential matches 
the size of the matter potential over an extended period of time. The transition continues as long as the neutrinos 
can change flavor in such a way that the potential matching is possible. Once the neutrinos can no longer keep up 
with the matter potential, then the transformation ceases. 

Whether a MNR region occurs and whether it results in flavor transformation depends upon the configuration of 
the emission surfaces. We examine two qualitatively different types of self interaction potentials that can arise from 
mergers. We point out that while there exist two different self interaction potentials, which result in distinct outcomes, 
the oscillation phenomena are both well described as MNR transitions. Our study is presented in this paper as follows: 
First we discuss representative disk configurations in section m then we discuss oscillation calculations in section Ell 
We present neutrino oscillation calculation results in section 113 and the results of nucleosynthesis calculations in 
section 13 We then conclude in section ED 


II. ACCRETION DISK CONFIGURATIONS 

Similar to core collapse supernovae, the mergers of compact objects like two neutron stars or a black hole and a 
neutron star release vast amounts of energy in the form of neutrinos and antineutrinos. However, unlike supernovae, 
these systems begin with a composition of almost entirely neutrons. After the collision, the material is heated and the 
electron fraction increases. As a consequence, the resulting emission of antineutrinos is _greater than the emission of 
neutrinos. Electron neutrinos and antineutrinos are thought to dominate the flux [D-H, [ill . [38l| and other flavors may 
not be trapped at all [s^ . Simulations of compact object mergers show small amounts of mu and tau neutrinos with 
a luminosity of up to 1.5 x 10®^ erg/s [ij. This luminosity is about 1/10**^ of the luminosity of electron antineutrinos. 

Guided by these results, we generate a set of representative disk models and calculate the neutrino oscillation 
pattern above them. The disk models used in the following calculations are inspired by those in [T^ liol - ld^ , for a 
merger of a 2.5Mq black hole and a 1.6Mq neutron star that forms a black hole with mass 3.85Mq and spin parameter 
a = 0.6. For ease of calculation, the disks are taken to be geometrically thin. We use two general types of models: 
one with all three types of neutrino emission coming from a disk with a single maximum radius, and a second where 
each the three types comes from a unique disk with its own radius. 

Each of the three kinds of neutrinos (electron neutrinos, electron antineutrinos, and all other flavors of neutrinos 
and antineutrinos) are taken to have different temperatures. These temperatures and radii are listed in Table |D for 
our single radius model and Table HD for our multiple radius model. We further vary the amount of Vr, and Vt 
within each type of model. 

We follow a test neutrino as it leaves each disk model at 45° from the plane of the disk. The test neutrino starts 


a 

Ta (MeV) 

Ra 

(cm) 


6.4 

4.5 

X 

10® 


7.4 

4.5 

X 

10® 

V<L 

7.1 

4.5 

X 

10® 


7.4 

4.5 

X 

10® 


TABLE I: Single Radius Model: Parameters for neutrino emission from a single surface. Fluxes are taken to be 
thermal Fermi Dirac fluxes for electron neutrino and electron antineutrinos. For mu and tau flavor neutrinos and 
antineutrinos, the flux is taken to be the thermal flux rescaled by an overall parameter specified in each example. 


a 

Ta (MeV) 

Ra (cm) 

Ve 

5.9 

5.2 X 10® 


9.9 

variable 

Ue 

7.8 

3.9 X 10® 

^fi,T 

9.9 

variable 


TABLE II: Multiple Radius Model: Parameters for neutrino emission from several surfaces. Eluxes are taken to 
be thermal Eermi Dirac fluxes for electron neutrino and electron antineutrinos. For mu and tau flavor neutrinos the 
fluxes are taken to be the thermal fluxes as well, but the radius of the emission is taken to be different as specified in 

each example. 
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FIG. 1: Trajectory of the test neutrino, f = {x,y = 0, z). 

above the disk and follows a radial trajectory outward, which we show in Fig. [T] The trajectory follows the mass 
outflow as it leaves from the disk. We use a parameterized outflow velocity, u, as in [i^. The velocity depends on 
the acceleration of the material, (3, and the eventual velocity of the material at infinity, Vcx, 

Rinner \ , 

1 -: ( 1 ) 

where R is the distance from the center of the disk and Rinner is the initial position of the material, which we take 
to be Rinner = 2.0 X 10® Cm. As the neutrino travels along the trajectory, the neutrino will feel potentials based on 
its position. The matter potential, due to coherent forward scattering with electrons and positrons, is computed from 
the number density based on the outflow model, which assumes a constant mass outflow rate. The outflow model, 
with s/k = 50, P = 2.0 and Voo = 0.1c, yields net electron number density, Ne{t), where t parameterizes the position 
along the trajectory. The position t = 0 is the start of the trajectory. 


III. CALCULATIONS 


We calculate the flavor transformation of neutrinos and antineutrinos as they travel along the trajectory which we 
show in Fig. [I] The evolution of neutrinos and antineutrinos is computed through the S matrices, which in the flavor 
basis are governed by. 


t-S{E) = (HviE) + Heit) + S{E). 

at 


The vacuum Hamiltonian is given in the flavor basis by, 

f-A2iiE) 0 

0 A2i(£;) 


Hv{E) — U2'3{02'i)Ul3{9i3)Ui2{9i2) 


\ 


0 

0 I ul2.(.9i2)ul^{6ii>,)ul^{922,), 

(A3i(£’) + A32{E)) 


( 2 ) 


(3) 


where Aij{E) = (rrif — m^)/(4£') and the Uij(9ij)s are the unitary matrices that take the Hamiltonian from the 
mass basis where E[v{E) is diagonal to the flavor basis. The matter potential influencing neutrino oscillations is 
Ve{t) = \/2GFNe{t), where Gp is the Fermi constant. The matter potential results in Hamiltonian contribution, 

He{t), 


He{t) 


( Ve{t) 
0 

V 0 


0 0 
0 0 
0 0 


(4) 
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Neutrinos also feel a potential from interacting with other neutrinos often called the self interaction potential. This 
potential is computed similarly, but because the neutrinos are not isotropically distributed above the disk, the test 
neutrino will feel a non trivial geometric effect from the other neutrinos. The resulting Hamiltonian contribution is 







S*{t,E)p{t,E)S^{t,E)) dE. 


(5) 


Initially S is the identity matrix, but evolves according to Eq. [5] as the neutrinos oscillate. The matrices, p and p 
take the form. 


/ ^.AE)c.At) 0 0 

p{t,E)=V2GF 0 KiE)C,At) 0 

V 0 0 cj>,AE)c,At) 

(^vAE)c,At) 0 0 

p{t,E)=V2GF 0 KiE)G^At) 0 

V 0 0 ^,AE)c,At) 


( 6 ) 


where Ga{t) is the geometric contribution, and the (f>aiE) is the flux for the a = disk. The geometric 

contribution was derived in [l8l| for a disk in general. Since we are calculating trajectories like those in Fig. [1] which 
make a 45° angle with the plane of the disk, we can reduce this general expression to 


Ca{t) 



G(r, x)rdr, 


where 


C{r,x) 


ttAx^ 

\/2 



(7) 

( 8 ) 


with Ra is the radius of the disk, E is the (anti)neutrino energy, I = {x — rA + x^, to = (x + r)^ + x^, £{M) is the 
complete elliptic integral of the second kind, and x is the distance from the center along the plane of the disk. 

We take the fluxes of neutrinos and antineutrinos to have the spectrum of the Fermi Dirac flux with zero chemical 
potential so that for electron neutrinos or electron antineutrinos with temperature, T, the flux is 




gc 


E^ 


27r^(?ic)^ 1 + ’ 


(9) 


where 5 = 1 is the spin parameter. For Vr{i’T) with temperature, T, the flux is also taken to have a Fermi 

Dirac spectrum with zero chemical potential. However, for the case of the single radius model, the flux is rescaled by 

fo, 


ME) 


. gc E^ 


( 10 ) 


where a = or W- For the multiple radius model, the flux is not rescaled so /o = 1 and the amount of zz^, 

zz^. Ft and Ft flux is entirely determined by the disk radius and temperature. Initially, the system has the same flux 
of zz_/i. Ft, Ffx and Ft from each flavor and these neutrinos all have the same emission geometry. Thus the density 
matrices can be written as. 


p{t, E) =V2 Gf 


p(t,E) =V2 Gf 


( M{E)CMt) 

0 

V 0 

( M{E)CMt) 

0 

V 0 


0 

M(E)CMt) 

0 

0 

M(E)CMt) 

0 


0 

0 

MiE)CMt) 

0 

0 

M(E)CMt) 


( 11 ) 
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At t = 0, the part of the Hamiltonian that comes from the neutrino-self interaction, (equation [S]) , has only one 
non-zero element which is the ee**' element. We define this to be W(t), 

K(i) =v.M-v,At) 

=V2Gf r (C., {E) - (E)) dE. 

do 

Sometimes it is instructive to examine the unoscillated potential, i.e. what the self-interaction potential would be if 
no oscillations were to occur. We will use V^, for this purpose. Of course, as neutrinos oscillate, the self interaction 
potential does evolve. The S matrices of Eq. [S]gain off diagonal components and therefore so does H,^^. The potential 
in the component deviates from Eq. 1121 and we define 


Eo.c = {H,,)ee - Tr{H,,)lZ. 


(13) 


Before calculating the neutrino flavor transformation above our two types of disk configurations, we compute the 
unoscillated self interaction potential, 14 (t), using the values for Ra and Ta from Tables HI and HIl We show the results 
in Figs yia\ and [5al for the magnitudes of all the potentials, A 12 , IA 32 I, Veit), and |14(t)|. Fig. [2a]is computed for 
the single radius model (Table |T| and Fig. [Ha] is computed for the multiple radius model (Table |lT]) . The vacuum 
potentials, A 12 , and IA 32 I, are plotted with a dark green dotted line and light green finely dotted line respectively. 
These vacuum potentials in the figure have a neutrino energy of F = 15 MeV. The self interaction potential, |14(t)| 
is plotted in a dashed light blue line. The horizontal axis in all plots is progress along the neutrino trajectory in cm. 
Since we use trajectories that are at a 45 degree angle the progress along the trajectory is ^ \/2x. 

Close to the disk, at the start of the trajectory, the geometric effect on the self interaction potential varies little. 
At these distances, where the scale of the disk is large compared to the distance of the neutrino from the disk, 
the geometric contributions of neutrinos and antineutrinos are similar: Cp,, ~ and roughly constant. The hotter 
antineutrino disk will contribute more flux ((/)p^ > ) and this higher antineutrino emission causes the self interaction 

potential to start out negative. Far from the disk, the geometric contribution to the self interaction potential has the 
form 


Gait) 


HK - Rfnner) 

256V2x(t)4 


(14) 


The turnover between these two behaviors is determined by the disk radius, Ra- In the single disk model of Fig. [5^ 
this turnover happens at about Xtum ^ 10® cm because all flavors have the same geometry, with Ra = 4.5 x 10® cm 
for all a. 

In the multiple disk model of Fig. I3al a similar turnover in the potential would occur at about Xtum ^ 10® cm as 
well. However, because the disks for electron neutrinos and antineutrinos have different radii, there is an additional 
feature. At Xsym ^ 10® cm, the neutrino self interaction potential for the multiple radius disk shows a symmetric 
point where the self interaction potential goes through zero, at the place where the electron neutrino and electron 
antineutrino contributions to the self-interaction potential are the same. This symmetric point is a consequence of 
the nontrivial geometric contribution to the self interaction potential (Eq. [T]) and the differing electron neutrino and 
electron antineutrino radii and temperatures. 

Again, the electron antineutrino disk is hotter than the electron neutrino disk, and the local emission surface 
produces more flux. 


/ <^p„(F)dA< / ^,AE)dE, (15) 

do do 

regardless of the local behavior of the geometric factors, C'p^(t) and C'p^(t). Close to the disk surface, the disk sizes 
are less important in determining the self interaction potential (Eq. 1121) than the temperatures since the geometric 
contributions are essentially the same and constant. Therefore the self interaction potential will be negative and stay 
roughly constant close to the emission surface as can be seen in Figs. l^andlHal 

While the flux of the electron antineutrinos always dominates over the flux of the electron neutrinos everywhere 
on the trajectory, farther from the disk surface, as the disk sizes become important, the self interaction geometric 
contribution associated with the neutrinos becomes larger than the geometric contribution associated with the an¬ 
tineutrinos, Gjj^it) < C'p^(t). This occurs because the radii, Ra, that set the scale of the potentials’ turnover are 
different: Rp^ < Rp^, so the geometric contribution for the electron antineutrinos begins to decrease sooner than the 
geometric contribution for the electron neutrinos (Eq. [HI) . The symmetric point at about Xsym ^ 10® cm is the 
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place where the two potentials are the same magnitude and sum to zero. Before this point is negative and after it 
becomes positive. 

The matter potential, which is shown as the light blue line in Figs. and [3^ depends on the mass density of 
material in the outflow and the electron fraction. The matter close to the disk where the outflow trajectory would 
be vertical stays at a relatively constant density, which we take to be close to 10^° g/cm^. As the outflow enters the 
radial trajectory, where we begin our calculation, the matter expands, dropping the density to roughly 10 ® g/cm®. 
After this turnover, as the mass streams far from the disk, it continues a decline in density, calculated assuming a 
constant mass outflow rate as described in Sec. El 

Before embarking on the neutrino flavor transformation calculation, we can use the (unoscillated) potentials plotted 
in Figs. [2a| and [3a] to identify regions where various types of oscillation physics may take place. In the single disk 
model. Fig. [3^ at the beginning of the trajectory, the matter potential, in a dark blue line, and the size of the self 
interaction potential, in the light blue line, are large. The large self interaction potential corresponds to synchronized 
neutrino oscillations where all modes behave roughly the same way, even though they have different energies. After 
the synchronized region, the matter potential and the neutrino self interaction potential have the same magnitude at 
2 X 10^ cm creating a matter neutrino resonance region. At about 10® cm where the neutrino self interaction scale 
and the vacuum scale are roughly the same size, there is a nutation region, and at about 3 x 10® cm, an MSW regime. 
The multiple disk model, Fig. [3^ has these regions of interest at more or less the same places and an additional 
matter neutrino resonance region at Xsym ^ 10 ® cm due to the symmetric point. 


IV. RESULTS 


We now wish to consider the oscillation of neutrinos for the models discussed in the previous two sections. We 
perform the calculations outlined in section uni for each of our models with varying contributions from Vt, 
and Vt- We take the vacuum mixing parameters to be consistent with the current values of the Particle Data Group 
HI, 6 »i 2 = 34.4°, 013 = 9°, 023 = 45°, A 12 A = 7.59 x 10"® eV^ and IA 23 AI = 2.43 x 10"® eV^. We use the inverted 
hierarchy although the matter-neutrino resonances are nearly hierarchy independent and we find very similar results 
for these transitions in the normal hierarchy. The self interaction Hamiltonian couples the neutrinos and antineutrinos 
of different energies together and our integration keeps track of 800 different neutrino energies between 1 and 101 
MeV. The result of the calculation is the S matrices, which we use to find the flux weighted survival probabilities. 




{P) = 


J^<P,AE)dE 


(16) 


The energy dependent survival probabilities and transition probabilities come from the S matrix. 



— 

(E)-- 

= \See 

iE)[ 



(E)-- 

= \S^.e 

m\ 

P.. 

— 

{E)-- 

= \Sre 

m 

P^e 

— 

(E)-- 

ICO 

II 

{E)\ 

P^. 

— 

{E)-- 


■XE)\ 

PP. 

— 

(E)-- 


■XE)\ 


(17) 


The evolution matrix, S, for the antineutrinos obeys an equation similar to Eq. dH). 

The capture rates for electron neutrinos and antineutrinos on neutrinos and protons respectively are approximated 
as 


Tr cGp 
^ _ TTcGp 


(4 + 3ci) ^ (r, E)P,^^,^ [E) + (r, E)P,^^,^ {E) + (r, E)P,^^,^ (A)) {E + Qf dE 

(4 + 3ci) (r, E)P,^^,^ {E) + 4^ (r, E)P,^^,^ {E) + 4^ (r, E)P,^^,^ (A)) {E - Q)^ dE 


( 18 ) 
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where Q = nin — rup is the nucleon mass difference and Eq = (m.„ — rup + nie) is the threshold energy for electron 
antineutrino capture. The fluxes, 4>'g^{r,E), decrease appropriately as the distance from the disk increases (this 
geometric effect is described in @ and is different from (f>a{E)Ca)- In addition to plotting survival probabilities we 
will plot the ratio of the capture rates in the case with neutrino oscillations to the case without neutrino oscillations. 
In the absence of oscillation, the capture rates, A° 


and 

A? , are calculated 

i/g 7 



=1 

P^.- 

^uAE) 

=0 

Pu.- 


=0 

P^.- 


=1 

P^.- 


=0 

Pp.- 


=0, 


(19) 


for all energies, E. 


A. Single Disk Models 


We first present the results of the calculations described in Sec. m using the single disk model discussed in Sec. 
im We calculate the neutrino oscillation pattern resulting from these potentials and show the results in Figs I2b[ 
and [2d] for different amounts of emitted mu and tau neutrinos and antineutrinos. 

The potentials in Fig. [5^ shape the oscillation regions that the neutrinos enter by their relative sizes. The matter 
potential and the self interaction potential become the same size at 2 x 10^ cm. We identify this position as a 
standard MNR region. Standard MNR transitions have been explained in [l^ and has similarities with the neutrino- 
antineutrino transformation discussed in [d^. In a standard MNR transition, the self interaction potential from Eq. 
(US changes so that it matches the size of the matter potential throughout the transition: |t4| ^ |I4sc|- This means 
that the on-diagonal component of H stays near zero, i.e. 

+ {Hv)ee 0 (20) 


If the system is sufficiently adiabatic, a standard MNR transition will occur: neutrinos and antineutrino oscillate 
in just such a way that the potentials maintain the resonance. The oscillated Vosc changes to enforce Eq. [501 As 
discussed in the behavior of this resonance transition is well approximated by a single monoenergetic neutrino 
and antineutrino model. In this context, [T^ derived the diabaticity criteria that comes from matching the timescale 
of transition determined by the scale height, r = \d{Ve/V,^^)/dl\~^, of the matter and self interaction potentials, 

dll = rlog , (21) 


with the timescale determined by the capacity of system to change flavor. 


Sh 


a 

A sin 20Sext ’ 


where 


a = 


YeM 

VuM 


( 22 ) 


(23) 


and Sext is the average difference between the y components of Neutrino Elavor Isospin vectors for the neutrino and 
antineutrino , which must be less than (1 -|- a)/2. Below this value, 6 li = SI 2 as long as 9 is sufficiently large. For 
the single disk model, the scale height is Tsingie ^ 7 x 10® cm at the standard MNR region, where, a ^ 1.4. These 
values mean that in order for a transition to occur associated with the mixing angle must be 6* > 3 x 10“^. 

The measured mixing angle, 613 ^ 0.2, is well within that range. The standard MNR regions in our models should 
be sufficiently adiabatic to sustain transitions. 





Survival Probability Survival Probability Potential (erg) 



Position (cm) 
(a) 



10^ 10* lo’ 10* 10^ io‘° 
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10* 10* lo’ 10* io‘' io‘“ 


Position (cm) 

(b) 



10* 10* lo’ 10* io‘' io‘“ 


Position (cm) 

(c) 


Position (cm) 

(d) 


FIG. 2: Single Disk Model: All flavors of neutrinos and antineutrinos are emitted from a disk of the same size. 
The neutrino and antineutrino temperatures differ, as in Table HI The horizontal axis in all plots is progress along 
the neutrino trajectory in cm. Fig. \2ai Potentials entering Hamiltonian from electrons as in Eq. |4] and from 
neutrinos as in Eq. [12] in the absence of oscillation. All other plots: The top panel shows the flux weighted electron 
neutrino survival probability, (P) in red solid lines, and the flux weighted electron antineutrino survival probability, 
(P), in dashed amber line. In the bottom panel, we show the relative capture rates of the electron neutrinos and 
antineutrinos. The relative capture rate of the neutrinos is a ratio of the electron neutrino capture rate when 
oscillations are taken into account, to the electron neutrino capture rate when oscillations are not present, 
and is shown as the dark blue line. The relative electron antineutrino capture rate is the analogous ratio, Ap^/A?^, 
which we show in a light blue dashed line. Fig. [EB: No mu or tau neutrinos are emitted from the disk. Fig. [23: Mu 
and tau neutrino fluxes are rescaled; /o=0.35 relative to their blackbody fluxes. Fig. \2di: Mu and tau neutrino fluxes 

are rescaled; /o=0.75 


We now look at our numerical calculations to see if the transitions indeed occur. The results of these calculations 
are shown in Figs. I2b[ I2cl and l2dl If there are no mu or tau neutrinos and antineutrinos, the transition results in flux 
weighted survival probability for electron neutrinos dropping to nearly zero and the flux weighted survival probability 
of electron antineutrinos returning to nearly one, as in Fig. I2bl This is the characteristic behavior of a standard 
MNR described in M- The relative capture rates of electron neutrinos are shown in solid dark blue and electron 
antineutrinos are shown in light dashed blue lines in the lower panel. Because no mu or tau neutrinos are present 
initially, these relative capture rates track closely to the weighted survival probabilities. 

Increasing the initial amount of mu and tau neutrinos by a small amount has little effect on the survival probabilities, 
but a significant effect on the capture rates. This can be seen in Fig. |2c|where a small quantity of mu and tau neutrinos 
are now emitted from the disk. The number flux of each of these types of neutrinos and neutrinos is slightly under 
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40% of the electron antineutrino flux. We see from a comparison of the top panel of Fig. I2bl and Fig. [23 that the 
addition of the mu and tan neutrinos at this modest level changes the survival probabilities very little. However, 
a comparison of the bottom panels shows that the relative capture rates have changed. Both the relative capture 
rates of electron neutrinos and of the electron antineutrinos begin at 1, like the weighted survival probability. During 
the MNR, the capture rate of electron neutrinos drops only to about 70%. The small change in capture rate occurs 
despite the fact that the MNR results in a strong transition from initially electron flavor neutrinos to other flavors, 
because mu and tau neutrinos also transform at the MNR to electron flavor. 

Increasing the initial amount of mu and tau neutrinos by a larger amount than in Fig. |23 has a more dramatic 
effect. In Fig. I2d[ where the mu and tau neutrino flux is 85% of the electron antineutrino component, the top panel 
shows that neutrinos initially in electron flavor remain in electron flavor, even as they pass through the MNR region. 
The bottom panel of Fig. |2d] is consistent with this failure of the MNR to result in transition; the capture rate of 
electron neutrinos remains the same, just as if there were no MNR region. This behavior occurs generically when 
the mu and tau fluxes approach the same level as the electron neutrinos and antineutrinos. During MNR transitions, 
the neutrinos and antineutrinos transform in such a way so that the self interaction potential matches the matter 
potential, Eq. 1201 But since a neutrino flavor transformation involves a trading of electron type neutrinos with mu 
and tau type neutrinos, if there are similar numbers of electron neutrinos as other types then the self-interaction 
potential, Eq. 1131 cannot change much and no transformation occurs. 

After the MNR region, the system passes through a nutation region, which happens at about 10® cm in Figs. I2b[[^ 
and l2dl Nutation oscillations are a high neutrino density phenomenon that have been deeply studied in the supernova 
setting, e.g. [1^. They have also been studied in the disk setting, They occur when the self interaction 

potential approaches the size of the vacuum scale and a required initial condition is that the electron neutrinos and 
antineutrinos begin approximately in their flavor eigenstates. During the transition the neutrinos and antineutrinos 
are said to be “locked” in that their transition probabilities are linked. The transitions are largest in the inverted 
hierarchy which is the case for the calculations presented here. In Figs. 1111123 and Hd] the nutation region occurs 
after the matter neutrino resonance region which would be a typical expectation in a disk system. These figures 
show that nutation oscillations occur only in the case where MNR transitions do not, i.e Fig. Hd] This is because 
the MNR transition moves the electron neutrinos and antineutrinos out of their initial eigenstates, and the required 
initial conditions for a nutation oscillation are not met. We note that nutation oscillations have been shown to depend 
upon multiangle effects which are expected to be large in a disk setting. The authors of [13 discussed a formalism by 
which to use a single angle calculation to find nutation oscillation survival probabilities that are robust in a merger. 

A Mikeyev Smirnov and Wolfenstein (MSW) region occurs when the matter potential becomes the same size as 
the vacuum scale. In the models of Figs [2bll23 and I2d[ this region is at about 4 x 10® cm. At the densities where 
MSW regions occur, the self interaction potential tends to be small and the neutrinos and antineutrinos of different 
energies act almost independently. The energy dependent oscillations will then occur only if the system is sufficiently 
adiabatic. In all three examples, we see that the MSW regime results in transitions of the electron antineutrinos to 
other flavors. These transitions of electron antineutrinos would not be possible in the normal hierarchy where instead 
electron neutrinos would transform. 


B. Multiple Disk Model 

In the multiple disk model, neutrinos of different flavors are emitted from disks of different radii as given in Table 
im There are three different emission surfaces, one for electron neutrinos, one for electron antineutrinos and one for 
all other flavors. The potentials for this model are shown in Fig. jjaj and the different regions of oscillation behavior 
were discussed in Sec. imi In this section we discuss the results of the multiple disk model calculations, which are 
shown in Figs. |3bl [33 andlSdl 

Initially, the multiple disk models behave in the same was as the single disk models and have a region of synchronized 
behavior close to the disk, so that the survival probability is essentially one in the region before ^ 10® cm. However, 
after this point the multiple disk models start to exhibit a differences with the single disk models because they enter 
the symmetric MNR region which does not exist in the single disk model. In the symmetric MNR region Vosc from 
Eq. [T3| becomes the same size as 14 ■ As discussed in section uni in the multiple disk models, the difference in the 
changing geometric factors causes the self interaction potential 14, that would otherwise have large magnitude, to 
sweep through zero. We call this a symmetric resonance region because it occurs close to the point where a = 1. 
The transitions in this region are similar to those observed in [l8l | and the phenomenology of these transitions has 
similarities to the standard MNR M- In all three of Eigs. [30123 andlSdl both neutrinos and antineutrinos undergo 
a nearly complete transformation. 
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FIG. 3: Multiple Disk Model Different flavors of neutrinos and antineutrinos are emitted from disks of different 
sizes, as in Table HIl i.e. = 5.2 x 10® cm and = 3.9 x 10® cm. The horizontal axis in all plots is progress 
along the neutrino trajectory in cm. Fig. I,9fll Potentials entering Hamiltonian from electrons as in Eq. |4]and from 
neutrinos as in Eq. [T2]in the absence of oscillation. All other plots: The top panel shows the flux weighted electron 
neutrino survival probability, {P) in red solid lines, and the flux weighted electron antineutrino survival probability, 
(P), in dashed amber line. In the bottom panel, we show the relative capture rates of the electron neutrinos and 
antineutrinos. The relative capture rate of the neutrinos is a ratio of the electron neutrino capture rate when 
oscillations are taken into account, to the electron neutrino capture rate when oscillations are not present, A®^, 
in a dark blue line. The relative electron antineutrino capture rate is the analogous ratio, Ap^/A®^, which we show in 
a light blue dashed line. Fig. [fR - No mu and tan neutrinos or antineutrinos are emitted. Fig. [2E' Mu and tau 
neutrinos and antineutrinos are emitted from a disk with radius of 1.8 x 10® cm. Fia \3d[ Mu and tau neutrinos and 

antineutrinos are emitted from a disk with radius of 2.4 x 10® cm. 


The mu and tau neutrinos can prevent the MNR transition if their flux becomes roughly the same as the electron 
neutrinos and antineutrinos. This case was shown in the single disk model in Fig. I2dl However, in all the examples 
considered as multiple disk models, the mu and tau neutrinos are either non-existent as in Fig. I3bl or emitted from 
a smaller disk than the I'e and Pe, Ru^ = 1-8 x 10® cm in Fig. [5R and = 2.4 x 10® cm in Fig. I3dl Therefore 
they do not give roughly the same contribution as the electron type neutrinos and antineutrinos at the position of 
the resonance and the transition can proceed. 

While the presence of a relatively modest number of mu and tau neutrinos does not effect the survival probability 
of the electron neutrinos and antineutrinos at the symmetric matter neutrino resonance, they do alter the capture 
rates. This is again because the electron neutrinos don’t simply “oscillate away”, instead they exchange places with 
the mu/tau type. For example, the lower panel of Fig. I5R shows in the blue solid line the relative capture rate of 
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electron neutrinos, which drops from one to only about 0.6 at the symmetric matter neutrino resonance transition. In 
the lower panel of Fig. I3d[ where the mu and tan contributions are a little larger, the electron neutrino capture rate 
is higher than had there been no transformation at all. This is due to the higher energy of the mu and tau neutrinos 
at emission and the energy squared dependence of the neutrino capture cross section. 

After the symmetric matter neutrino resonance transition, the neutrinos encounter a standard MNR region at about 
2 X 10^ cm in Figs I3bl ISH and I3dl Whether or not a standard MNR transition occurs depends on the state of the 
system as it enters the standard MNR region. Transitions cannot occur unless the resonance condition, Eq. [501 is 
fulfilled and for this to happen, the neutrino self interaction potential must be negative. While the potential, W begins 
negative close to the disk, the changing geometric factors cause a relative sign change in the potential (see Eq. USD 
so that 14 is positive at the position of the standard MNR. However, 14 is what the self-interaction potential would 
be in the absence of oscillations. The system has already undergone a transformation at the symmetric MNR and 
this transition produces an additional change in the sign so that the actual self-interaction potential, I4sc, is negative. 
Thus if a symmetric MNR transition occurs, the system is set up up favorably for a standard MNR transition. Since 
in all our examples a symmetric MNR transition occurs, we expect a standard MNR transition as well. The change 
in survival probability at the standard MNR region can be seen in the top panels of Eigs. |3bl and [3dl to follow 
the typical behavior of a standard MNR transition. Similar to the single disk model, the standard MNR transitions 
have a non trivial effect on the capture rates as can be seen in the bottom panels of these figures. 

After the MNR regions, the neutrinos encounter the nutation region at about 10® cm. As in the single disk models, 
Figs. I2bl and [5H the system is no longer in the original flavor eigenstates and the nutation region results in no flavor 
transformation. The neutrinos then enter an MSW region after several times 10® cm, where some flavor transformation 
takes place in Figs [33 |3cl and [3dl Also, as in the single disk case, since we are using the inverted hierarchy, the 
electron antineutrinos in the multiple disk case transform to other flavors at the MSW region. 


V. NUCLEOSYNTHESIS 

The symmetric MNR region shown in Fig. occurs quite close to the disk and thus will impact the element 
synthesis in outflowing material from the inner disk regions. This material starts out as primarily free nucleons, 
where neutrons far outnumber protons, and as it moves outward the composition evolves via the weak interactions 

Vf. + n ^ p + e~ (24) 

De+P ^ n + e'^. (25) 

The approximate rates for the forward reactions are given in Eq. 1181 In the nucleosynthesis calculations described 
in this section we include also the weak magnetism contribution |46j to these rates. Since the disk emits more 
antineutrinos than neutrinos and the antineutrinos tend to be hotter, neutron-rich outflows are favored; the outflows 
may also retain some of the neutron-richness of the disk. As the material expands and cools, the protons are quickly 
bound into alphas, and of the weak reactions above only the top forward reaction, Vf, + n ^ p + e~, continues to 
operate. The protons thus produced are promptly bound into alphas as well. At this stage of the nucleosynthesis 
neutrinos act to increase the number of seed nuclei and reduce the number of free neutrons available for capture on 
the seeds. This is called the ‘alpha effect’ [13,14^, and it limits how far in A the nucleosynthesis can proceed. Most 
calculations of merger disk outflow nucleosynthesis [1, H, [l^, [H, IH) IHi I13| favor production of weak (80 < A < 120), 
rather than main {A > 120), r-process nuclei. The MNR can potentially alter this conclusion. 

The influence of an MNR transition close to the emission surface on outflow nucleosynthesis was first pointed out for 
a collapsar-type disk in . A collapsar disk emits primarily electron flavor neutrinos and antineutrinos, so when the 
MNR transition occurs it is (from the perspective of the nuclear matter in the vicinity) as if the neutrinos disappear. 
If this happens in the region of the outflow where alphas are forming, the alpha effect can be cut off completely. With 
fewer seeds formed and more free neutrons remaining, a vigorous main r process can result 

Exactly how the MNR will influence merger outflow nucleosynthesis depends on the relative amounts of electron 
and mu/tau flavors emitted. If there is little mu/tau emission, we expect an effect similar to that described above 
for the collapsar case. This is shown in the dark/light blue lines in Figs. |3|and|51 Here, we start with the neutrino 
emission from the multiple disk example, described in Sec. IIVBI and calculate the element synthesis as in [I^, [l^ 
along the outflow trajectory described in Sec. [Hi The top panel of Fig. 0] shows the final abundances for the case 
with no neutrino oscillations, compared to the solar r-process abundance pattern. Primarily A ^ 80 — 90 neutron-rich 
nuclei are produced. The material starts out very neutron-rich, as shown in the middle panel of Fig. [SJ but the alpha 
effect limits the resulting nucleosynthesis to a weak r process. When the neutrino oscillations illustrated in Fig. 3 are 
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FIG. 4: Nucleosynthesis resulting from the multiple disk models. The vertical axis shows abundance and the 
horizontal axis shows the mass number A of elements produced. Both Panels: The black pluses show scaled solar 
r-process residuals. Top Panel: Production from the multiple disk model with no oscillations is shown in black. 
Bottom Panel: Production from the multiple disk model with oscillation that includes no initial mu and tau 
neutrinos is shown in dark blue. Oscillation calculations where mu and tau neutrinos have a flux relative the 
electron antineutrinos at 5% (light blue), at 10% (green), at 20% (yellow) and at 65% (red) are shown. 


included in the calculation, we find the MNR can radically change this picture. Fig. I3bl shows that the symmetric 
MNR region occurs at a point about 10 km along the outflow trajectory. If there are little to no mu/tau neutrinos 
to oscillate into the electron flavors, the electron neutrino capture rate will drop steeply in the MNR region. This 
is shown for example cases with little (light blue lines) to no (dark blue lines) mu/tau emission in the top panel of 
Fig. O For the outflow conditions considered here, MNR occurs just before the alphas start forming, as depicted in 
the bottom panel of Fig. [5j Thus, the alpha effect is completely removed by the MNR —fewer alphas form (bottom 
panel of Fig. El), more free neutrons remain (middle panel of Fig. [5]), and a robust main r process results (bottom 
panel of Fig. SI). 

If there is appreciable mu/tau emission from the disk, then the electron capture rates will not drop so steeply in the 
MNR region and may even increase, as described in Sec. IIVI and shown in Figs. IScland lSdl To explore this effect on 
the element synthesis we repeat the oscillation and outflow calculations described above for increasing percentages. 
Since these are multiple disk calculations, we adjust the radius of the emission surface of the mu/tau neutrinos so 
that it corresponds to 10%, 20%, and 65% of mu/tau emission as compared with electron antineutrino emission. The 
results are shown in Figs. |4] and [5] by the green, yellow, and red lines, respectively. The symmetric MNR influences 
both the electron neutrino and antineutrino fluxes, effectively swapping them with hotter but weaker mu/tau fluxes. 
Thus for a short time before alphas form, the balance in the forward reaction rates of Eq. [25] is adjusted by the 
MNR and the balance of protons and neutrons correspondingly shifts as shown in the middle panel of Fig. [5] In the 
10 % case, the neutrino capture rates still decrease, but there are enough electron neutrinos after the transition that 
a modest alpha effect occurs. This means that there are not enough neutrons for the r-process to get all the way to 
the highest mass number nuclei, although the nucleosynthesis does proceed beyond the A ^ 80 r-process peak region 
and some A ^ 130 peak nuclei are produced. In the 20% case, the neutrino capture rates are relatively unchanged 
but the antineutrino capture rates are reduced, so that the overall composition is more proton rich than in the no 
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FIG. 5: Top Panel: Neutrino capture rates from the multiple disk models as a function of position along the outflow 
trajectory. Middle Panel: Neutron mass fraction as a function of position along the outflow trajectory. Bottom 
Panel: Alpha particle mass fraction as a function of position along the trajectory. All Panels The multi disk model 
with no oscillations is shown in black. The multiple disk model with oscillation that includes no initial mu and tau 
neutrinos is shown in dark blue. Oscillation calculations where mu and tau neutrinos are included at 5% are shown 
in light blue, at 10% in green, at 20% in yellow and at 65% in red. 


oscillation case. This combined with a robust alpha effect results in a sharp reduction in A ~ 80 nucleosynthesis 
compared to the no oscillation case. With 65% mu/tau emission, the MNR results in faster neutrino capture rates, 
and the composition shifts to roughly equal numbers of neutrons and protons before alphas begin to form. Here only 
iron peak nuclei result. 


VI. CONCLUSIONS 

Compact object mergers emit large numbers of neutrinos with a flux that is initially composed of more electron 
antineutrinos than neutrinos. This makes this environment a prime candidate for matter neutrino resonance transi¬ 
tions. There are two types of relevant neutrino emission configurations: one where all flavors of neutrinos are emitted 
from the same surface, and one where neutrinos are emitted from different size surfaces. The latter is more in line 
with the results of recent compact object merger simulations. We find that two types of matter neutrino resonances 
occur. Both configurations show a standard MNR, which produces a transition where electron neutrinos wind up 
transformed into other flavors and electron antineutrinos wind up back in their initial eigenstates. We also find a 
symmetric MNR, where both electron neutrinos and antineutrinos are nearly completely transformed. The type of 
transition is a consequence of the behavior of the neutrino self-interaction potential which depends sensitively on the 
balance of r'e and DgS- 

The size of the initial contribution of mu and tau type neutrinos to the flux, as compared with electron antineutrinos 
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and neutrinos, is crucial to determining whether a transition occurs or not. Comparable contributions will shut off 
the transition, but such large and fluxes are not currently predicted. 

Future calculations of neutrino transport in mergers will be instrumental in determining what oscillation pattern we 
can expect above real world mergers, both because it is the balance of and i>eS that determines the initial potential, 
and also because the mu and tau fluxes determine the flexibility available to the system for a MNR transition. 

Matter-neutrino resonance transitions have a significant impact on wind nucleosynthesis, because they occur when 
the neutrino self-interaction potential and matter potential are approximately balanced. If such transitions occur, 
they often occur close to the emission surface of the neutrinos, where the neutrinos still have large enough flux to 
affect the neutron to proton ratio, and nuclei in the outflow have not yet begun to form. 
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